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ABSTRACT 



Context. Strong gravitational lensing is a unique tool that can be used to model with great accuracy the inner mass distribution of 
massive galaxy clusters. In particular, clusters with large Einstein radii provide a wealth of multiply imaged systems in the cluster 
core. Measuring the redshift of these multiple images provide strong constraints to determine precisely the shape of the central dark 
matter profile. 

Aims. This paper presents a spectroscopic survey of strongly lensed galaxies in the massive cluster lens Abell 1703, displaying a 
large Einstein radius (28" at z = 2.8) and numerous strongly-lensed systems including a central nVzg-like configuration. 
Methods. We used the LRIS spectrograph on Keck to target multiple images and lensed galaxy candidates, and use the measured 
spectroscopic redshifts to constrain the mass distribution of the cluster using a parametric model. 

Results. The spectroscopic data enable us to measure accurate redshifts for 7 sources at z > 2 , all of which in good agreement with 
their photometric redshifts. We update the identification of multiply imaged systems by discovering 3 new systems and identifying 
a radial counter image. We also report the discovery of a remarkably bright ~ 3.6 L* /-band dropout at z = 5.827 in our mask that 
is only moderately magnified by the cluster (// ~ 3.0 + 0.08). The improved parametric mass model, including 16 multiple systems 
with 10 spectroscopic redshifts, further constrain the smooth cluster-scale mass distribution with a generalized NEW profile of best-fit 
logarithmic slope a = 0.92 + 0.04, concentration C200 = 4.72 + 0.40 and scale radius = 476 + 45 kpc. The overall RMS in the image 
plane is 1.3". 

Conclusions. Using our strong-lensing model, we predict a large scale shear signal that is consistent with weak-lensing measurements 
inferred from Subaru data out to 4 Mpc h~^. Together with the fact that the strong-lensing modeling requires a single dark matter 
clump, this suggests that Abell 1703 is be a relaxed, unimodal cluster. This unique cluster could to be probed further using deep 
X-ray, SZ and dynamics analysis, for a detailed study of the physics in a relaxed cluster. 
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1. Introduction 

Massive galaxy clusters with large Einstein radii (Re ^ 2C) of- 
fer a wealth of information about their central mass distribution, 
by displaying large numbers of multiple images in their strong 
lensing regions. The most striking case i s the cluster Abel l 1689, 
with Re ~ 50^^ (lBroadhurst et al.ll2005l iLimousin et"aD r2QQ7). 
which holds the largest number of strong lensing constraints in 
a single field. This gives a great opportunity to map inner mass 
profiles to a high level of precisi on. Parametric cluster-lens mass 
models have been widely used (iKneib et alJ 119961: ISmith et aP 



moaels nave been widely 

iRichard et al.ll2QQ7l) . but to construct accurate mass mod- 
els and obtain their fiducial parameters, accurate redshift infor- 
mation is necessary. Indeed, one redshift of a set of multiple im- 
ages enables us to estimate the absolute total mass within the 
Einstein radius. Adding many more redshifts allow to measure 
the full central mass profile accurately. Even though photomet- 
ric redshifts offer an acceptable level of precision for such mod- 



eling (iBroadhurst et al.ll2QQ5l: iBradac et all '2008V spectroscopic 
redshifts are the most reliable way to confirm the identifications 
and find new multiply-imaged systems. Lowering down the red- 
shift error on the multiple systems with spectroscopy has an- 
other great interest, as one can use strong-lensing to constrain 
the cosmological parameters w and based on purely geomet- 
rical considerations (Link & Pierce 1998; Golse et al. 2002). 
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On the theoretical side, IBroadhurst & Barkanal (l2008l) re- 
cently compared the measured lensing parameters for a sam- 
ple of clusters with large Re, with predictions from the standard 
ACDM cosmology, and report a 4-cr discrepancy in the pres- 
ence of such large Re for the corresponding virial mass. Their 
result, based on a sample of massive haloes derived from the 
Millennium simulation, includes a full treatment of the projec- 
tion eflTects, which largely bias this sample towards alignments of 
the triaxial dark-matter haloes along the line of s ight. This bias 
is also pointed out by Oguri & BlandfordI (l2008h . who investi- 
gated a full-sky Monte-Carlo realizations of massive clusters. 
They show that although the population of clusters with large 
Einstein radii is scarce, it is possible to find Einstein radii as large 



as 5C in the whole sky, under the assu mptions of a WMAP5 
cosmology. The main diff erence bewe en TBroadhurst & Barkanal 
(I2QQ81) and .Oguri &Blandfo^(l2QQ 8') is the latter' s treatment of 
the sample as a whole- sky survey, which is far more representa- 
tive of the observed population than the Millennium simulation, 
which only covers a small volume of -5000 Mpc^.This releases 
any disagreement between the presence of large Re and A-CDM. 

Among the clusters with large Re.^Q have started to study 
in lLimousin et al ] (I2008h . hereafter LOS, the mass distribution of 
Abell 1703, one of the richest clusters discovered in the Sloan 
Digital Sky Survey (SPSS) and showing strong gravitational 
lensing (iHennawi et al ] l2008h . In this previous work, we have 
identified 13 systems forming highly-magnified multiple im- 
ages, including a central ring. Based on the spectroscopic red- 
shift of this central feature and photometric redshifts for the 
remaining images, all the l ensing systems have been well re- 
produced by a sing le NFW (lNavarroetal.1119971) profile for the 
dark matter, making Abell 1703 relatively simpler compared to 
bimodal lensing clusters (e.g. Eliasdottir et al. 2007 for Abell 
2218, Richard et al. 2007 for Abell 68, Verdugo et al. 2007 for 
MS2053). We therefore decided to use it as a more reliable probe 
for comparing the results of strong and weak-l ensing measure- 
ments in clusters with large Einstein radii from lBroadhurst et al.l 
(I2OO8). 

We present in this paper results from a spectroscopic sur- 
vey conducted at Keck using the LRIS instrument, aiming at 
confirming the photometric redshifts, multiply images identifi- 
cation, and precisely constraining the mass distribution in Abell 
1703. We obtained new spectroscopic redshift measurements for 
most of the systems identified by L08 (10 out of 13), in good 
agreement with their corresponding photometric redshifts. We 
also report of an exceptionally bright /-band dropout at z ~ 6, 
forming a magnified single image close to the region of multi- 
ple images. We further identify three new systems of multiple 
images and constrain the lensing model, as well as compare the 
weak-lensing extrapolation of our model to the recently mea- 
sured weak-lensing profile of Broadhurst et al. ( 2008). 

We present in Section 2 the observations and data reduction 
of the spectroscopic data. Section 3 gives the redshift measure- 
ments, spectra and the main results from the mass modeling. We 
discuss these results in Section 4. Throughout the paper, we used 
a standard A-CDM model with = 0.3, Da = 0.7 and = 70 
km s"^ Mpc "\ whenever necessary. In this cosmology, V on 
the sky is equivalent to a physical distance of 4.244 kpc at the 
redshift z = 0.28 of the cluster. All magnitudes are quoted in the 
AB system. 

2. Observations 

2.1. Imaging and photometric redshifts 

The current study is based on the same imaging data described 
in L08, based on Hubble Space Telescope and Subaru observa- 
tions of Abell 1703 covering broad-band filters from B to H. 
These images were previously use d to derive photometric red- 
shifts with the software Hyperz (Bolzone lla et al.l 12000). We 
briefly summarize the characteristic parameters, central wave- 
length and depth, of these images in Table [T] 

2.2. Keck spectroscopy 

We ha ve used the Low Resolution Imager and Spectrograph 
(LRIS. IOke et al.|[T995l) on the Keck telescope to perform multi- 
slit observations of the Abell 1703 cluster field. We designed 



a multi-slit mask containing 32 slits of 1.0'' width to include 
as many of the multiple systems previously identified in L08 as 
possible, with a few tilted slits following the geometry of long 
arcs (Fig.[T]). In the case of fainter or less reliable identifications, 
independent images of the same system were observed in sep- 
arate slits. Additional slits in the mask included new multiple 
systems candidates, cluster members and background galaxies 
as estimated through their photometric redshifts. 

The LRIS instrument was set up with a 6800 A dichroic iso- 
lating the blue arm, with a 300 lines/mm grism blazed at 5000 A, 
from the red arm, with a 600 lines/mm grating blazed at 7500 A. 
This ensures a high throughput over the wavelength range 3500 
to 9500 A and good spectral resolution in the red to resolve sky- 
lines as well as [Oil] emission line doublets in the redshift range 
0.85 - 1.40. This instrument setup is also well-matched with the 
photometric redshift estimates of the targetted sources. We ob- 
tained 2 hours of exposure time on 2008 May 10 in photometric 
conditions and very good seeing (0.7"). 

2.3. Data reduction and redshift measurements 

Data reductio n of the spectr a was performed using the Python 
version of the iKelsonI (l2003l) reduction scripts, which off'er the 
advantage of processing the images in their distorted framework. 
This helps to reduce noise correlations, in particular for the case 
of tilted slits. We performed standard reduction steps for bias 
removal, flat-field correction, wavelength calibration, sky sub- 
traction and cosmic-ray rejection, and used observations of the 
standard star BD-h33-2642 obtained on the same night to derive 
the flux calibration. 

Close-ups on the targets and extracted spectra are presented 
in Fig. [2] and [3] We measured the redshift for 33 sources in the 
entire LRIS mask. The spectra show either Lyman-of (in absorp- 
tion or emission) with additional ultraviolet absorption lines in 
the blue part, or a resolved doublet of [Oil]) in the red part 
of the spectrum. In two multiple systems (10/11 and 16), we 
observed these images in multiple slits: redshift measurements 
were obtained by identifying significant absorption lines after 
stacking the relevant exposures. For each spectrum, the average 
redshift value is obtained from the peaks of the main spectral 
features identified, while the corresponding error is taken from 
the spectral dispersion. Additional uncertainties generated by the 
accuracy of the relative and absolute wavelength calibrations, of 
about 1.1 and 1.5 A, respectively, were quadratically added to 
yield the final redshift errors. A confidence class Zdass, ranging 
from 1 to 4, was assigned to each measurement according to the 
prescription of iLe Fevre et all (|l995): this corresponds to a prob- 
ability level for a correct identification of 50%, 75%, 95%, and 
100%, respectively. A specific value of 9 is used when only a 
single secure spectral feature is seen in emission. 

Table [2] summarizes the redshift measurements Zspeci confi- 
dence class Zciass and the spectral features used to derive the red- 
shift for all objects detected in the LRIS mask, including those 
additional targets serendipitously falling into the slits. For the 
majority of sources located within the ACS field of view, we 
report in the same table the corresponding photometric redshift 
Zphot- We find a good agreement between individual Zspec and 
Zphot values, with a typical error < \zphot - Zsped >= 0.15 or 
< \^ph£il^\ >=0.047, and no catastrophic measurements. This 
Strong consistency is mainly owing to the use of space-based im- 
ages of very similar PSF size, together with ground based near- 
infrared images obtained with MOIRCS in excellent seeing con- 
ditions (OA''). 



3. Results 



3.2. Spectroscopy of a remarkably bright i-dropout 



3.1. Multiple images redshifts and identifications 

Figure [T] presents a color image of the central region of Abell 
1703, combining the B, V and Z ACS filters, where we marked 
in white the location of the 1 1 multiply imaged sources identified 
in our previous work (LOS). Our spectroscopic data enabled to 
measure the redshifts for the following 8 systems (displayed as 
white squares in Fig.[T]): 

- System 1 at z = 0.8889 forms a very highly magnified ring- 
shaped configuration with 4 images close to the BCG and 
a less-magnified 5th image to the south-west. Although its 
redshift was previously derived in L08, we obtained a new 
spectrum with the same mask in much better conditions. 

- System 3 at z = 3.277 is located to the north of the 
BCG, with 3 images: 2 bright merging images and a fainter 
counter-image to the west. 

- Systems 4/5 at z = 1.9082 are two individual regions identi- 
fied in the same source which forms a cusp configuration of 
3 images. 

- System 6 at z = 2.360 is a similar configuration as Systems 
4/5, but closer to the center of the cluster. 

- System 7atz = 2.983 is again a cusp configuration like 
Systems 4/5, but further from the center of the cluster. 

- Systems 10/11 at z = 2.627 are two individual regions iden- 
tified in the giant tangential arc at the south-east of the BCG. 
A fainter counter-image is located to the west. The new mass 
model presented here enable us to associate the radial feature 
reported in L08, very close to the BCG (next to the ring con- 
figuration of System 1) to be a 4th image (10.4) associated 
with this arc. 

- System 15 at z = 2.355 is a typical Einstein- cross configura- 
tion of 4 images surrounding the BCG. We note that systems 
6 (z = 2.360) and 15 (z = 2.355) have similar redshifts and 
are therefore very close to one another in the source plane, 
with a projected separation of 67+6 kpc. 

- System 16 at z = 2.810 is another Einstein- cross configura- 
tion located at a slightly higher redshift, but found very close 
to System 15 in the image plane due to projection eff'ects. 

Additional sources 20 and 21, which were included in the 
mask as a possible new multiple system, show a very similar 
spectroscopic redshift z = 1.279 with [Oil] doublet of emission 
lines. However, they are located too far from the center of the 
cluster to be multiply imaged at this redshift. 

Apart from system 1 , all the other systems are located in the 
redshift range 1.9 < z < 3.3. Using the new spectroscopic red- 
shifts into the mass model (see Sect. l3.3l) , we have computed the 
predicted region of multiple images for a redshift z = 6 (delim- 
ited with a dashed line in Fig.[T]). As we expect all the multiple 
images present in the HST optical bands to lie within this re- 
gion, we use it to update our systematic search for multiple im- 
ages based on the photometric redshifts catalog. We find three 
additional systems complementing the list presented in L08. We 
display them as red circles in Fig. [T] and report the corresponding 
positions and photometric redshifts in Table O 

- System 12 is again a typical cusp configuration of 3 images 
located to the north, and "following" the nearby Systems 4/5, 
7, 8 and 9. 

- System 13 and 14 are two further Einstein-cross configura- 
tions of 4 images surrounding the BCG. System 14 is located 
near the Systems 15 and 16 in the image plane, while System 
13 is closer to the cluster center. 



One of the sources in the LRIS mask, identified as 23 in the 
spectroscopic catalog (Table [2|, deserves particular attention. 
Located to the south-west of the cluster center (green triangle 
in Fig.[T]), it has been selected for its very red i - z = 2.9 + 0.3 
color and bright z = 24.30 + 0.07, while being undetected blue- 
ward of the / band in the other deep ACS images. The NICMOS 
andMOIRCS photometry / = 23.93+0.07 and if = 23.75+0.10 
further constrain its spectral energy distribution to be compati- 
ble with a photometric redshift Zphot = 6.0 + 0.42 (Figure (H top 
and middle panels), and the / - z color make it a very strong i- 
d ropou t photometric candidate according to the iBouwens et al.l 
(l2006h selection criterion. 

The LRIS spectrum for this object covers the wavelength 
range 3800-9230 A, corresponding to 2.12 < z < 6.59 for 
Lyman-Qf. We detect a 5 cr- significant emission line at /I = 8300.5 
A (FigurelH lower panel), with an integrated flux / = 2.5 + 0.4 x 
10"^^ ergs s"^ cm"^, and an observed equivalent width W> 50 
A when assuming a 5(T upper limit on the underlying contin- 
uum.This line is slightly resolved in the spectral direction, with 
a possibly asymmetric profile. In order to estimate this asymme- 
try, we use the weighted skewness statistics 5"^, which is based 
on the third momentum of the flux distribution (seen as an 
array of size N) and has been used as a se lection criterion for 
Lyman-Qf emitters (LAEs. iKashikawa et al.ll2006) . The value of 
5'w, measured in Angstroms, is defined as: 

= Yj^Xi - xf (^10,. - Aio,b) (1) 

where (x/) and (fi) are the arrays of coordinates and fluxes, of 
size N, / = Yjf fi^ and (x, a) represent the mean and standard 
deviation of the array (x/), respectively. The weight factors /lio,r 
and /1 10,/? correspond to the wavelengths where the flux drops to 
10% of it s peak value on the red and blue sides of the line emis- 
sion (see IKashikawa et aDl2006l for a more detailed discussion 
about the use of 5'^). In our case, we measure /lio,r = 8305 A 
and /lio,^ = 8298 A on the spectrum of the dropout, and a corre- 
sponding 5"^ = 5. 6+1. 4 A. The error on 5"^ has been measured 
using a bootstrap resampling of the fluxes. The skewness value 
is slightly superior to the critical value 5^ = 3 A which has been 
us ed as a lower lini it for selecti ng single line emitters as LAEs, 
in lKashikawa et al.1 (^2006) and Shimasaku et al.l (l2006h . The sin- 
gle emission line is therefore likely to be Lyman-or at a redshift 
z = 5.827, quite compatible with the photometric redshift pre- 
diction. 

At his redshift, we compute a magnification factor and as- 
sociated error yu = 1.2 + 0.03 magnitudes, derived with the im- 
proved strong lensing model (see Sect. 13.3.21) . We note that al- 
though this object is a single image, it is located very close to 
the region of strong lensing, allowing an accurate estimate of the 
magnification with our model. 

After correcting for this magnification, the source is both 
magnified and intrinsically bright (Huniensed = 24.9), as there 
are only 9 of such sources at z ~ 6 brighter than z = 25.5 
AB found in the UDF, HUDF and GOODS fieds (iBouwens et al.1 
I2006h . This object is similar in redshift and intri nsic magnitude 
to the spectroscopically confirmed /-dropouts by IStanway et al.l 
(2004) at z = 5.78 and z = 5.83, or the slightly magnified 
object at z = 5.515 found in the cluster RDCS 1252.9-292 
(iDow-Hvgelund et al.l[2QQ5h . Based on the latest constraints of 
the luminosity function at z ~ 6 (IBouwens et al.| [2008). the un- 
lensed magnitude corresponds to an 3.6 L* galaxy. We also note 




Fig. 1. ACS color image of Abell 1703 (combination of the HST filters F450W, F606W and F850LP), showing the location of all 
the multiply-imaged systems included in this study. The white cross at the center of the image marks the location of the brightest 
cluster galaxy, which has been subtracted for clarity. The dashed line outlines the limit of the region where we expect multiple 
images up to z = 6. The systems presented in LOS are shown in white, and the new systems presented in this work are shown in red. 
Squares identify the multiple images with spectroscopic redshifts.The green triangle located to the lower right marks the location of 
the /-band dropout presented in Section [T2l 



that the very red i - z = 2.9 color i s much redder than for any 
of the 9 dropouts with z < 25.5 from iBou wens et all (l2006l) . de- 
tected in both the / and z bands, and where the reddest color was 
found to be / - z = 2.4. Thanks to the magnification and the 
combination of filters at this redshift, the / - z color allows to 
measure the average depression factor < D >=< 1 - /obs//int > 
at z ~ 6, between the observed and intrinsic fluxes shortward 
of Lyman-Qf, due to line blanketting. The corresponding value 
< Z) >~ 0.96 is a lower limit, due to possible contamination of 
the red side of the line in the / band filter. This value is close (but 
slightly superior) t o the predicted value at z = 6 (< D >= 0.92) 
when following thelMadaul(ll995l) prescription. 

Although the current photometry, probing only the rest- 
frame UV at z = 6, does not allow to derive physical param- 
eters on this source (such as age, reddening, stellar mass), we 
can measure the UV spectral slope ft, defined as /i oc /l^. Using 
the J - H = 0.18 color, we derive JS ~ -1.5, a value redder 



than a spectrum flat in fy (fi = -2.0), as can be seen in Fig. |4l 
Although this val ue is typical of Lym an-break galaxies at z ~ 3 
OS = -1.5 ± 0.4. lAdelberger & Stei del 2000 ), it is redder tha n 
for the sample of 27 /-dropouts studied by Sta nway et al.l (l2005l) . 

who derived JS 2.0 based on / and H p hotometry, or the sani - 

ple of lensed z-band dropouts presented in lRichard et al.l (l2008l) . 
The use of IRAC observations probing the longer wavelengths 
would allow us to understand whether this redder color is due to 
an old stellar population or aneffect of reddening, similar to the 
work done by 'Eg ami etaP (12005^) on a highly lensed galaxy at 
z ~ 6.8 discovered in AbeU 2218 (iKneib et al.ll2QQ4h . 



Fig. 3. Continuation of fig.|2] 
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3.3. Strong Lensing 
3.3.1. Modeling method 

Our starting point in modeling of the dark matter distribu- 
tion in Abell 1703 is the mass model presented in deta ils by 
LQ8, where we used the LENSTOOlHI public software ( Julio et al] 
l2Q07l) to constrain a parametric mass model with the identi- 
fied multiple systems. The model is optimized through the new 
Bayesian Markov chain Monte-Carlo (hereafter MCMC) sam- 
pler, described in detail in I Julio et aO (i2QQ7ii) . This process uses 
the observational constraints (positions of the multiply imaged 
systems) to optimize the parameters describing the mass distri- 
bution by matching the location of each image of a given system 
in the source plane. 

The contribution of the stellar mass from the central 
galaxy is parametrized by a Pseudo Isothermal Elliptical Mass 



Distribution(PIEMD) using the profile derived from the photom- 
etry. The PIEMD profile, whic h has been w idely used for the 
modeling of cluster- sc ale (S mith et al .1120051) and galaxy-scale 
(iNatarajan et al.|[T998h haloes, assumes each dark matter clump 
can be parametrized by a central position, ellipticity 43 , position 
angle 0, central velocity dispersion ctq, and two characteristic 
radii: a core and a cut radius. The total mass is proportional to 

2 

Using this parametrization, the PIEMD mass density S(x, y) 
takes the form: 



^0 



^cut 


r 1 


^cut ~ ^ core 





1 



(2) 



where x and y coordinates are oriented along the position 
angle 6, = [(x - x,)/(l + €)f + [(y - yc)/(l - €)]\ e = e/(2 - 



|http : //vivivi . oamp . f r/cosmology/lenstool/] 



^ We use here the definition e = 1 - Z?/^ of the ellipticity, where a and 
b are the semi-major and semi-minor axis of the ellipse, respectively. 
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Fig. 4. A bright /-dropout found in the field of Abell 1703. The 
top panel shows the individual detections in the ACS, NICMOS 
and MOIRCS image. The object is resolved and elongated along 
the shear direction, with an angular size 0.45''. Middle panels: 
best fit of the spectral energy distribution with the photometric 
redshift software (left) , and corresponding redshift probability 
distribution (right). Bottom panels: VIZ ACS color image with 
the location of the LRIS slit (left), and close-up on the extracted 
spectrum where we detected an emission line at 8301 A (right). 
The grey region highlights the location of a bright OH skyline. 



e), and (Xc, yc) is the center of the mass distribution (see also 
ISmithet al.l (,2005) for more details). 

The cluster- scale dark-matter compone nt is niodele d with a 
generalized NFW (hereafter gNFW) profile fchaolll996h includ- 
ing a central logarithmic slope a\ 



p(r) 



Pc^c 



(r/r,)-(l + (r/r,))(3--) 



(3) 



where is a scale radius, pc the critical density and Sc is 
related to the value of the concentration parameter C200 through 
the relation: 



200 



^200 



3 ln(l + C200) - C2oo/(l + C200) 



(4) 



Again, the spherical gNFW profile is generalized to an ellip- 
tical mass distribution using a relation similar to Eq. [2l A more 
detailed description of th e implementation of the gNFW profile 
in lenstool is provided bv lSand et"aD (l2008l) . 

Individual cluster galaxies were added as small-scale 
PIEMD perturbers based on their photometry and shape param- 
eters, using empirical scaling relations between their dynam- 
ical parameters (central velocity dispe rsion and sc a le rad ius) 
and their luminosity, assuming the Fab er & JacksonI (Il976h re- 
lation and a constant mass-to-light ratio for each galaxy (see 



ICovone et al .120061 for more details). The same parameters of one 
particular cluster galaxy to the north (852) were optimized indi- 
vidually to reproduce the configuration of the surrounding sys- 
tems 3, 7, 8, 9 and 12. 

We replaced the photometric redshifts kept as free param- 
eters in L08 by their spectroscopic equivalents for the 7 rele- 
vant sources (9 systems). The precision of the new model, which 
uses the same parametrization as L08, can be estimated using the 
RMS of the location of multiple images in the image plane, de- 
fined as : 




xpred 



)2 -h (yobs^. ^ ■ 



ypredy 



(5) 



where (xobs j^kjobs j^k) and (xpredj;y^,ypredj;y^) are the ob- 
served and predicted location of image j in system k, respec- 
tively. Using the exact same parametrization, we find a value of 
(T/ = 1.3''similar to the result obtained in L08 (1.4''). This pre- 
cision is typical of strong lensing works using a similar number 
of multiples images (Richard et al. 2007; Eliasdottir et al. 20Q7|), 
and smaller by a factor of ~ 2 from the one obtained for Abell 
1689 r2.87". lLimousin et al.ll2Qb71 

We summarize the best fit parameters of the new model in 
Table |4l The main gNFW component shows a concentration 
C200 = 4.72^Q 4^, a scale radius = 476.9^42 g ^P^' inner 
slope a = 0.92^Q Q4. The geometrical parameters (center, orien- 
tation and ellipticity) are ver y simi lar to the ones measured on 
the cD galaxy (see also Sect. 13.3.41 below). In comparison with 
the results from L08, the strongest variations in the best fit pa- 
rameters are found in the gNFW profile, which appears to be 
more concentrated (larger C200 and smaller scale radius r^), but 
both results are compatible at the 3cr level. The addition of the 
new constraints (spectroscopic redshifts, new multiple systems) 
while keeping the same model parametrization, did not reduce 
the uncertainty in the best fit parameters, even if both models 
reproduce the constraints on the multiple images with a simi- 
lar precision. This argues for a remaining degeneracy within the 
gNFW parameters, which cannot be completely disentangled by 
a mass model purely based on strong lensing. 



3.3.2. Redshift predictions and magnifications 

One of the main benefits of the Bayesian approach is that the 
MCMC optimization provides a large number of models which 
sample the posterior probability density function of all the pa- 
rameters ( Julio et al. 2007). We can use these different realiza- 
tions to estimate the average value and associated error for any 
given parameter or combination of parameters. In our case, we 
estimate the redshift and associated error for each of the new 
multiple systems 12, 13 and 14. These values are reported in 
Table [51 We also computed the magnification factors yu and as- 
sociated error for every source in our LRIS mask (Table ^ as 
well as the other multiple systems (Table [3]) . The magnification 
factors range from 0.01-0.05 mags for single images to about 4.8 
magnitudes for the source 1.1. 

3.3.3. Weak-lensing predictions 

Our mass model, solely based on strong lensing constraints out 
to 50" (the region of multiple images), can be used to predict the 
tangential shear induced on background galaxies up to larger 
scales. The main observable with weak-lensing is the reduced 
shear g = y/(l - k) where y is the gravitational shear and k the 
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Fig. 5. Tangential reduced shear g+ predicted with the strong 
lensing model based on constraints in the central region (black 
error bar on top). The mean value and 3cr error is presented as 
a function of radius in the grey-filled region. We overplot the 
datapoints obtained by Broadhurst et al. ( 2008) in their weak- 
lensing measurements with Subaru. Extrapolation of our strong 
lensing model forms a very good fit to these observations upto 
r = 1000 arcsec, or 4 Mpc h~^. 

convergence. In order to compare our predictions with weak- 
lensing observat i ons of Abell 1703 with Subaru described by 
iBroadhurst et al.l (l2008l) . we adopt the definion of the tangen- 
tial reduced shear g+ from'Umetsu & B roadhurstl(l2008l) . i.e. the 
projection of the shear perpendicularly to the direction of the 
cluster center (see also lMedezinski et al]l2007l) . 

Using the MCMC realizations, we determine predictions for 
and its associated S-cr error as a function of angular distance r 
from the center of the mass distribution, assuming a distribu- 
tion of background sources at a fixed redshift Zs = 1.2 corre- 
sponding the the me an of the redshift distribu tion of sources in 
the observations from lBroadhurst et aP (12008 ^ Figure[5]presents 
these estimations against their weak-lensing measurements over 
the same range of r. We find that our strong-lensing model rep- 
resents a very good fit to the weak-lensing datapoints up to 4 
Mpc/z"\ with a of 4. 1/7 = 0. 58, even after extrapolating our 
mass model to a radial distance 20x larger than the region of 
strong lensing constraints. This argues again for the overall sim- 
plicity of the cluster, with a single NFW profile providing a good 
fit both for the strong-lensing and weak-lensing measurements. 

3.3.4. Ellipticity of the mass distribution and Einstein radius 

From the best fit geometrical parameters found during the opti- 
mization of the mass model (Table I?]), we notice a good agree- 
ment between the values for the ellipticity and position angle of 
the large scale gNFW profile (^nfw = 0.23, ^nfw = 63.9 degrees) 
and the same parameters of the cD galaxy (eco = 0.19, 6cd = 52 
degrees) as measured with SExtractor u .This is not really sur- 
prising, as this correlation has previously been reported both in 
strong-lensing analysis (e.g. Gavazzi et al. 2003 in MS2137), X- 
ray observations (Hashimoto et al. 2008) as well as numerical 
works (Dubinski 1998; Faltenbacher et al. 2005^ . 

More generally, the total mass distribution of dark matter is 
the sum of this gNFW component and the various PIEMD pro- 
files associated with the cluster galaxies, and the global elliptic- 
ity ^mass may have a diff'erent value. We use the convergence map 
K (which shares an identical geometry with the mass distribution) 
to measure this ellipticity emass posteriori by fitting elliptical 

^ The value of e^^pw given in LOS had an error in its definition and we 
update it to the correct one. 



contours with the IRAF routine ellipse. We keep the center 
of the ellipses fixed at the peak of the mass map and let the el- 
lipticity emass{(^) and position angle 6{a) vary as free parameters 
with the semi major axis a. In the central region (a < 15''), we 
find a best fit value emass ~ 0-23 (Fig. [51 in blue) dominated by 
the geometry of the cD and the gNFW profile. At larger distances 
emass^ci) increases to reach a maximum of 0.4 at a ~ 40''(Fig.[6l 
in green), due to an alignment of bright galaxies (as seen in pro- 
jection) to the north-west and the south-east, therefore along a 
similar direction as ^nfw- These galaxies contribute to increasing 
the value of emass compared ioe^pw^ this eff'ect being significant 
on the mass distribution up to the edge of the ACS field of view. 
The same ellipticity is seen on t he we ak-lensing reconstructions 
performed by Broadhur st et al.l (l2008l) . 

We also compute the effective Einstein Radius Re, defined 
as the radius r from the center of the cluster (located at the 
peak of the mass dist ribution) at which K(r < Re) = 1 
dBroadhurst & Barkanal 1200 8). We measure the value and error 
on Re for the redshift z = 2.627 of the giant tangential arc 
(formed with the systems 10 and 11), using the diff'erent MCMC 
realizations. We find Re = 28.0 ± 0.25 arcsec, a value smaller 
than one would measure based on the simple distance of the gi- 
ant arc (Fig.O. This value is close but smaller than the estimate 
of 32'' given by Broadhurst & Barkana (2008), who mentions 
an agreement of Re with the distance of the arc. This appar- 
ent discrepancy is due to the ellipticity of the mass distribution, 
as the giant arc is oriented along the major axis of the ellipti- 
cal mass distribution. This ellipticity is now well-constrained by 
the large number of spectroscopically confirmed multiple im- 
ages presented here, making the new measurement of Re more 
robust. 

We can use the same cumulative p robability distribution 
P(> Re) as Broadhurst & Barkanal (12008!) to test the agreement 
of this new Einstein radius with the ACDM Millennium sim- 
ulation. This probability is computed from the distribution of 
concentrations at a given virial mass, accounting for the vari- 
ous biases arising from selecting a population of lensing clus- 
ters, as well as the projection eff'ects. We find a higher prob- 
ability P=20% (instead of 7.9%) of agreement, making Abell 
1703 less discrepant than the other mentioned clusters (Abell 
1689, CL0024 and RXJ1347) for its virial mass. Together with 
the results from these 3 other clusters (P=8.5%, P=3.9% and 
P=13% respectively) the combined probability for these 4 clus- 
ters is 9e~^, or a 3.7cr discrepancy. However, as mentioned ear- 
lier, there is a very small number of such clusters observed in 
the whole sky in comparison with the size of the Millennium 
Simulation, and therefore it seems more plausible to reconcile 
the observations with whole-sky simulations. 

4. Abell 1703: a relaxed, unimodal cluster ? 

Our strong lensing analysis using a simple NFW component for 
the large scale dark matter distribution, after individual treatment 
of the central cD galaxy and small scale mass distributions as- 
sociated with the galaxy substructure in the cluster, is able to re- 
produce the large number of multiple images identified in Abell 
1703. Even if the overall RMS in the predicted location of the 
multiple images (1.3") is still quite larger than the precision in 
astrometric measurements with HST (typically 0.1"), we find no 
systematics as a function of the location of the multiples, that 
would argue for the need of secondary da rk matter clumps. This 
was the case in th e clusters Abe ll 2218 (iEliasdottir et al.ll2QQ7l) 
and Abell 68 (Ric hard et al.ll2007) which show a strong evidence 
for bimodality. 




Fig. 6. Convergence map distribution k for a source at z = 2.627, 
the redshift of the giant tangential arc (systems 10 and 11, lo- 
cated in red). By averaging this map inside increasing aper- 
tures, we measure the corresponding effective Einstein radius 
Re, slightly shorter than the distance to the giant arc due to the 
overall ellipticity of the mass distribution. The blue and green 
ellipses show the best-fit ellipses found on the overall mass dis- 
tribution, with an increasing ellipticity when moving from small 
to large radii (see Sect. 13.3.41 for details). 



Furthermore, adding new spectroscopic measurements on 7 
multiple sources has removed an identical number of free pa- 
rameters (i.e., the redshifts of the sources) from the analysis 
presented in LOS. However, we found no significant increase in 
the overall RMS of the multiple images compared to this previ- 
ous work, under the same assumptions. In addition, the strongest 
variations in the best fit parameters compared with the previous 
values presented in LOS are found for the gNFW profile, but both 
results are compatible within the 3 cr confidence level. The re- 
maining scatter in the positions of multiple images is very likely 
due to small-scale substructure in the dark matter distribution, 
which is not included here apart from the visible cluster galax- 
ies. One way to increase the complexity of this model would 
be to estimate small-scale deviations from the smooth mass dis- 
tribution, using adaptive-grid based models (Julio et al. 2009, 
submitted). 

Finally, the best fit parameters of our strong-lensing analysis 
have shown to be consistent with weak-lensing measurements. 
Our best fit concentration parameter C200 = 4.72^^4^, corre- 
sponding to Cyir = 5.S2^Q^^ with the second definition based on 
the virial radius rv,>, is somewhat lower than the value derived by 
iBroadhurst et al.l(l2Q0 S). c,„v = 9.92+^-^^, but marginally compat- 
ible at 2 Act level. Despite this smaller concentration parameter, 
the predicted shear profile provides an overall good fit to their 
observations: the main diff'erence is our smaller Einstein radius 
(which s ets the location of the first datapoint in Fig. \5\ as was 
done by IBroadhurst et all I200S). allowing for a smaller slope 
of the shear profile. We also note that the two shear measure- 
ments at large radii (where the signal gets weaker) are slightly 
lower than the predictions (but consistent within their error bars). 
Together, these diff'erences aff'ect the best fit concentration pa- 
rameter. The new Einstein radius measurement should be more 



reliable, since it is based on a large number of multiple images at 
the center. There seems to be a smaller discrepancy between our 
new strong tensing NEW parameters and the w eak tensing mea- 
surements, as it was the ca se for Abell 16S9 (IBroadhurst et al.l 
120051: iLimousin et al.l[2007h . 

Overall, we argue that these results provide further evidence 
for Abell 1703 to be a relaxed cluster, as was already suggested 
in LOS. Follow-up X-ray observations of the cluster core should 
allow to confirm again this assumption, and also to separate the 
X-ray gas component from the dark matter in the mass distribu- 
tionwhendoing the str ong-lensing analysis, following the work 
bv lBradac etaiT(l200Sl) in RXJ1347. Another independent mea- 
surement of the mass profile would be to study the dynamics of 
cluster galaxies up to large radii, giving another estimate of the 
virial radius and virial mass. Finally, we expect to find a larger 
sample of similar relaxed clusters from the ongoing LoCuSS 
survey (IZhang et al.ll200Sl) of clusters with lower masses (as in- 
ferred from their X-ray luminosity) at z ~ 0.2. By pursuing the 
same strong tensing analysis on this broader sample, we will 
probe a wider distribution of Einstein radii, ellipticities and con- 
centration parameters with less bias compared to available nu- 
merical simulations. 
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Table 1. Summary of the broad-band images used for the pho- 
tometry. The depth is given as the 5cr detection limit for a point 
source. 
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ID 



a (J2000.0) 6 (J2000.0) 



3.35+0-052 

3.35: 
2 25+ 



F775W jjL (mags) Zdass Features 



I. 1 
3.1 
3.3 

4.2/5.2 
6.1 
7.3 
10.1 

II. 3 
15.1 
15.3 

16.1/16.3 



198.77725 
198.76696 
198.75824 
198.76075 
198.77984 
198.75869 
198.78708 
198.76242 
198.76284 
198.78825 
198.78851 



51.81934 
51.83205 
51.82982 
51.82487 
51.82640 
51.82814 
51.81424 
51.80954 
51.81246 
51.82173 
51.82096 



0.8889+0.0003 
3.277 +0.002 
3.277 +0.002 
1.9082 +0.004 
2.360 +0.002 
2.983 +0.002 
2.627 +0.004 
2.627 +0.004 
2.355 +0.002 
2.355 +0.002 
2.810 +0.004 



2.59 
3.20 
3.10 
2.70 
2.44 
2.60 
2.75 



22.19+0.01 4.78+0.22 4 [Oil], [Nelll], H(5, Hy 

22.69+0.02 3.40+0.47 4 Lja(em), SilV, Sill 

25.34+0.05 1.66+0.22 4 Lj^Cem), SilV, Sill 

23.72+0.02 2.54+0.03 3 CIV, Sill, Znll 

23.56+0.02 1.96+0.04 4 Lj«(em), SilV, CIV, Sill 

26.81+0.10 1.87+0.02 9 Ly^i^m) 

21.40+0.01 2.69+0.65 3 Lj^Cabs), SilV 

23.55+0.02 1.73+0.02 3 Lj^Cabs), SilV, Sill 

24.89+0.05 2.26+0.10 4 Lj^Cem), CIV 

25.44+0.05 1.56+0.05 4 Lj^Cem), SilV, CIV 

25.52+0.06 1.61+0.02 3 Lj^Cabs), Fell 



20 

21 

23 

gal2 

gal4 

gain -a 

gall2 

gall3-b 

gall6 

gain 

gall8-b 
gall9-a 
gall9-b 



198.76659 
198.76000 
198.75612 
198.78811 
198.76596 
198.76772 
198.77770 
198.77518 
198.79701 
198.79065 
198.78816 
198.77100 
198.77063 



51.83626 
51.83451 
51.80736 
51.84276 
51.85902 
51.79678 
51.79462 
51.78871 
51.77606 
51.77097 
51.77364 
51.76841 
51.76935 



1.279 + 0.0005 
1.279 + 0.0005 
5.827 + 0.0015 
0.5632 + 0.0004 
3.2269 + 0.002 
2.2864 +0.002 
2.290 +0.002 
1.035 +0.001 
0.3906 +0.0005 
1.236 +0.0005 
0.5860 +0.0005 
0.2909 +0.0005 
0.3622 +0.0005 



1.10 

6.0+^-42 



•0.38 

550+0-^2^ 
3 44+0.09 

2 56+^'^^^ 

-0.252 

2 32"*" 

1 05+^'^^^ 



21.86+0.01 1.51+0.03 4 [Oil] 

22.57+0.02 1.75+0.02 4 [Oil] 

27.20+0.30 1.20+0.03 9 Ly^i^m) 

20.49+0.01 0.29+0.01 4 [Oil], K,H, Mgl, Mgll 

24.16+0.03 0.40+0.02 3 Lj^(abs), Sill, SilV 

22.86+0.01 0.85+0.02 4 Lj^Cem), Sill, CIV 

25.32+0.08 0.83+0.03 4 Lj^Cem), CIV, AlIII 

23.45+0.03 0.48+0.01 4 [Oil] 

0.10+0.01 4 K,H 

0.27+0.02 4 [Oil] 

0.19+0.01 4 [Oil], Mgll, K, H 

0.01+0.01 4 K,H 

0.06+0.01 4 [Oil], K,H 



gall 

gal3 

galll-b 

gall3-a 

gall4 

gall8-a 

gal20 

gal21 

gal22-a 

gal22-b 

gal23 

gal24 



198.78155 
198.77735 
198.76774 
198.77540 
198.78653 
198.78753 
198.79633 
198.78704 
198.77246 
198.77236 
198.75441 
198.75110 



51.83971 
51.85023 
51.79744 
51.78808 
51.78574 
51.77299 
51.76602 
51.76265 
51.75876 
51.76022 
51.86258 
51.86508 



0.1432 +0.0005 
0.087 +0.0005 
0.280 +0.001 
0.1762 +0.0005 
0.1011 +0.0003 
0.2705 +0.0005 
0.2766 +0.0005 
0.2781 +0.0005 
0.2410 +0.0005 
0.2766 +0.0005 
0.1765 +0.0004 
0.2683 +0.0005 



4 [Oil], H/?, [OIII] 

4 [Oil], K, H, Hyg, [OIII] 

4 [Oil], H/?, [OIII] 

4 K,H 

4 [Oil], H/?, [OIII] 

4 K,H 

4 K,H 

4 K,H 

4 [Oil], K,H 

4 K,H 

4 [Oil], K, H, H/?, [OIII] 

4 K,H 



Table 2. Overview of the sources included in the LRIS mask. From left to right: identification, astrometric position, spectroscopic 
redshift, photometric redshift estimate, ACS-F775W band photometry when available, magnification in magnitudes, redshift quality 
class (see Sect. 12.31) . main spectroscopic features. The mentions (abs) and (em) refer to a Ly^, line in absorption or in emission, 
respectively. The first part of the table corresponds to multiply-imaged system, the second part to background sources, the third part 
to cluster members and foreground sources. 
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(^JZUUU.UJ 


Zmodel 


Zphot 


r / / J W 


fi (mags) 


New multiple images 


10.4 


198.77173 


51.81906 


Zspec = 2.627 


~ 


24.77+0.08 


1.55+0.17 


12.1 


198.77456 


51.829946 


3.31+0.16 


3 33"^" 


26.68+0.10 


1.93+0.07 


12.2 


198.76568 


51.828658 




3 35+^-^^ 
-Q-21 


24.47+0.04 


1.87+0.16 


12.3 


198.75867 


51.826482 




3 16+0-24 


26.33+0.08 


2.17+0.04 


13.1 


198.78238 


51.821204 


0.839+0.007 


0.68!|| 

^ ^ ^ _i_n Aq 
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26.93+0.11 


1.48+0.03 


13.2 


198.76796 


51.818374 




24.18+0.04 


2.42+0.33 


13.3 


198.77204 


51.814973 




26.29+0.08 


2.16+0.19 


13.4 


198.76692 


51.815130 




25.12+0.05 


2.69+0.06 


14.1 


198.78662 


51.820542 


1.58+0.02 
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25.55+0.11 


1.63+0.02 


14.2 


198.76868 


51.820191 






25.40+0.17 


1.26+0.48 


14.3 


198.76447 


51.812192 






26.33+0.16 


2.13+0.03 


14.4 


198.77540 


51.811945 






26.30+0.15 


0.81+0.10 


Other multiple systems from LOS 


2.2 


198.77156 


51.81174 




2 31+0.46 

2 23+^"^^ 


25.29+0.08 


4.45+0.18 


2.3 


198.76970 


51.81186 


26.15+0.10 


4.37+0.16 


8.1 


198.77250 


51.83045 


9 -3 +0.004 
^•^■^■^-0.077 




24.95+0.06 


2.46+0.18 


8.2 


198.76608 


51.82949 


2 yy+O.!! 

9 79+8:20 
^•'^-0.12 


23.06+0.02 


2.19+0.10 


8.3 


198.75863 


51.82740 




25.07+0.06 


2.00+0.03 


9.1 


198.77176 


51.83030 


a 971 +0.031 
•^•^^ ^-0.124 




26.68+0.15 


3.09+0.20 


9.2 


198.76690 


51.82957 


9 QQ5+0-195 

;rQ-378 

3 00+0-37 


25.32+0.08 


2.48+0.19 


9.3 


198.75813 


51.82708 




26.92+0.12 


1.98+0.02 


Table 3, 


. Location and redshift estimate of three new multiply-image systems, as well as the radial feature (10.4) associated with 



the giant arc (Systems 10 and 11). From left to right: identification, astrometric position, redshift estimate from the lensing model, 
photometric redshift estimate, magnification in magnitudes. We also report the same values for the three remaining systems from 
LOS (Systems 2, 8 and 9), used in the mass model but lacking spectroscopic redshifts. The redshift of 10.4 has a spectroscopic 
measurement and therefore is not a free parameter in the model. System 14 and image 9.1 are too faint to provide a reliable 
photometric redshift estimate. 



Clump (x) (y) e 6 r (kpc) a C200 ctq (kms~0 Zmodd 

NFW -0.79+0.08 0.91+0.08 0.230+0.006 64.0+0.3 ^IG.^l^'l 0.92^0-05 4772^^ - = 

cD [0.0] [0.0] [0.19] [52] [25] 355.7:;o-3 

Galaxy 852 [19.0] [54.0] [0.11] [65.5] 98.7+1.8 320.9+3.5 

L, elliptical galaxy ^^-^-^il 202+^ 

Table 4. Best fit parameters used to model the mass distribution. The positions (x) and (y) are in arcsecs relative to the central 
galaxy (R.A.=13:15:05.276, Decl. =+51:49:02. 85) and oriented in WCS (north is up, east is left). Orientations 6 are measured in 
degrees and ellipticities (e) of the potentials are given as 1 - b/a, where a and b are the semi-major and semi-minor axis of the 
ellipse, respectively. Values in brackets were not optimized by the model. The radius (r) refers to the scale radius in the case of 
the gNFW profile, and the rcut radius in the case of the PIEMD profile (see LOS for details). 



